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Abstract 

We discuss the Schwinger mechanism in scalar QED and derive the 
multiplicity distribution of particles created under an external electric 
field using the LSZ reduction formula. Assuming that the electric field 
is spatially homogeneous, we find that the particles of different momenta 
are produced independently, and that the multiplicity distribution in one 
mode follows a Bose-Einstein distribution. We confirm the consistency 
of our results with an intuitive derivation by means of the Bogoliubov 
transformation on creation and annihilation operators. Finally we revisit 
a known solvable example of time-dependent electric fields to present exact 
and explicit expressions for demonstration. 
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1 Introduction 



A classic example of a non-perturbative tunneling phenomenon in quantum field 
theory is the decay of an electric field due to pair creation. The quantum vac- 
uum is full of virtual particle-antiparticle pairs (i.e. the Dirac sea), which can 
occasionally gain enough energy from the external field to become real. The 
decay (or persistence) rate of the QED (quantum electrodynamics) vacuum 
in the presence of an external electric field was first deduced from the imagi- 
nary part of the Heisenberg-Euler Lagrangian [1] and formulated in Schwinger's 
classic paper [2]. The phenomenon is commonly referred to as the Schwinger 
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mechanism (see rcf. [3] for a comprehensive review.) In the case of QED the 
coupling constant e is very small, and it is difficult in practice to achieve large 
enough electric fields; the probability of producing an electron-positron pair is, 
up to a prefactor, ~ exp[— nm^/ (eE)\, and thus a very strong electric field, 
E <~ m\/e ~ 1.3 x 10 18 V/m, is necessary to observe the phenomenon. To 
the best of our knowledge, the Schwinger mechanism in QED remains to be 
unambiguously observed experimentally. 

Pair creation from an electric field became a phenomcnologically much more 
relevant subject with the realization that the strong nuclear force is described 
by a gauge theory called QCD (quantum chromo-dynamics). A popular phe- 
nomcnological view of QCD with confinement is a description in terms of a 
chromoelectric flux tube connecting the color charges of the quarks. If these 
quarks are then pulled apart by their momenta, the string formed by the chro- 
moelectric field can decay via the Schwinger mechanism leading to the decay of 
the system into qq or color neutral mesons as a result of hadronization. In this 
case, the decay probability is characterized by ~ exp[— nm^/a], where the QCD 
string tension a ~ 1 GeV/fm is an energy stored in the chromoelectric flux tube 
per unit length. Applications of the particle production by the Schwinger mech- 
anism range from e + e~ annihilation [4, 5] to early models of relativistic heavy 
ion collisions [6, 7, 8, 9, 10, 11, 12]. Recent extensive studies on thermal hadron 
production as a possible manifestation of the Hawking-Unruh effect, that is an 
equivalent formulation to the Schwinger mechanism in curved space-time, is 
found in refs. [13, 14, 15]. 

The QCD coupling constant g, although asymptotically small, is not as small 
as the QED one at phenomcnologically interesting energies. Even more impor- 
tant is that the nonlinear dynamics of the gauge fields naturally leads, in some 
circumstances, to gauge fields that are parametrically large in the coupling, 

~ l/g. A prime example of such a situation is caused by the large occupa- 
tion numbers of gluonic states in high energy scattering. The transverse gluon 
density ~ Qs provides a typical energy scale Q s in such a system. There the non- 
linear interactions among bremsstrahlung gluons with small Bjorken's x lead to 
gluon saturation, which is most conveniently described as a coherent color field 
radiated by static (in light cone time) sources. This description is referred to as 
the Color Glass Condensate (CGC) (for reviews, see [16, 17, 18, 19]). The colli- 
sion of two objects whose wavefunction is characterized by Q s in the CGC for- 
malism achieves a field configuration of longitudinal chromoelectric and chromo- 
magnetic fields whose strength is also given by Q s . This transient state of matter 
containing strong longitudinal fields is known as the glasma [20, 21, 22]. In the 
case of QCD it is of course difficult to achieve the canonical model case of a con- 
stant electrical field. Generically, a WKB-type non-perturbative evaluation such 
as in refs. [13, 14] could be expected to be valid in a case where the field strength 
gA^ ~ Q s is much larger than the typical (inverse) time and spatial scales of the 
fields. A perturbative calculation, on the other hand, is also feasible in the case 
when g is small enough. The particle (gluons and quarks) production associated 
with strong CGC fields has been formulated based on the Lehmann-Symanzik- 
Zimmermann (LSZ) reduction formula [23, 24, 25, 26, 27, 28, 29, 30, 31] as well 
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as on the canonical formalism [32]. One of the aims of this paper is to establish 
a link, by a concrete example, between the formalism of LSZ reduction formulas, 
which is usually associated with perturbation theory only, and non-perturbative 
tunneling phenomena of the Schwinger mechanism. More concretely, we want to 
show that the LSZ perturbative framework automatically includes the particles 
produced by the Schwinger mechanism, provided the external field is properly 
resummed. Thus, this contribution does not need to be added separately by 
hand. 

For applying the mechanism of pair creation from a classical field to phe- 
nomenology one needs, in addition to the vacuum decay rate or spectrum of 
pairs, the whole probability distribution of multiparticle production. In many 
practical applications of the Schwinger mechanism, there has been a confusion of 
terminology, with both the formulas and the concepts of the vacuum decay rate 
(or persistence probability calculated by Schwinger) and the pair production 
rate. The difference between the two was recently nicely discussed in ref. [33], 
where it is interpreted as a result of temporal correlations between the produced 
pairs. In fact these two were clearly distinguished already in a classical paper 
by Nikishov [34] . In the case of the typical QED discussion, the pair production 
rate is extremely small, in which case the probability distribution of produced 
pairs cannot be distinguished from a Poisson distribution. This seems to have 
been the assumption used, without any further justification, also in many QCD 
phcnomenological applications (see e.g. [8] where this is very explicit). As we 
shall show explicitly in the following, this assumption is not true when the pair 
production is not strongly suppressed, as can typically be the case in e.g. heavy- 
ion collisions. Instead, the probability distribution of the produced pairs turns 
out naturally to be the appropriate (Bose-Einstcin or Fermi-Dirac) quantum 
one 1 . With explicit expressions for the probabilities to produce one, two, etc. 
particles, the distinction between the vacuum decay rate (related to the proba- 
bility to produce no pairs) and the pair production rate (the expectation value 
of the number of pairs produced), becomes obvious. The fact that there is a 
quantum statistical (BE or FD) correlation has long ago been realized by some 
authors as a requirement that should in principle be built into Monte Carlo 
event generators [37, 38, 39] . Also, the full computation of the vacuum decay 
rate should encompass all the multiparticle production processes, -because of 
unitarity-, including the quantum statistics. To our knowledge, however, an 
explicit derivation of how the BE or FD correlations arise from the Schwinger 
mechanism has been lacking. Besides, the multiparticle distribution has scarcely 
drawn attention in the context of the Schwinger mechanism, probably because 
of the absence of experimental access. In fact the spectrum of multiparticle 
production is quite informative and precise data of charged hadron multiplicity 
fluctuations are already available in p-p [40] and heavy-ion [41] collision experi- 
ments, where a negative binomial distribution gives a beautiful fit. Of course, to 
account for the high-energy experimental data, a simple treatment of spatially 

1 Note that the "inversion of spin statistics" discussed in [35, 36] refers to a formal expression 
of the vacuum decay rate as an integral over the BE or FD distribution function and not the 
actual probability distribution of produced particles. 



3 



homogeneous (i.e. constant in space) background fields is inadequate, and re- 
cently, it has been shown that an inhomogencous configuration forming a certain 
number of the glasma flux tubes leads to the negative binomial distribution [42] . 
This is beyond the scope of our current paper. 

In the following, we shall study the case of scalar QED in a time-dependent 
but spatially homogeneous external gauge field. We shall first introduce the 
model and derive the probability distribution of produced pairs using the LSZ 
reduction formula in sec. 2. Then, in sec. 3, we shall rederive the same results 
using canonical quantization and interpret the calculation of sec. 2 in terms 
of a Bogoliubov transformation. The discussion in sec. 3 to a large degree 
follows that of Tanji [43], but takes the additional mathematically simple step 
of actually writing down the whole probability distribution (see also [44] ) . Then, 
in sec. 4 we shall demonstrate how this procedure works in practice using an 
exactly solvable example (see also [45]) of a time dependent external potential, 
from which we can take both the constant field and short pulse limits. We note 
that, in all our discussions, we will solve the problem for a given external field 
without taking into account the interplay between the field and the produced 
particles which screen the external field, leading to plasma oscillation behavior 
in time [43, 46, 47]. 

2 LSZ derivation 

We here calculate the Schwinger mechanism in terms of the LSZ reduction 
formulas. To this aim we develop a slightly modified version of the Schwinger- 
Keldysh formalism to compute the generating functional of the particle and 
antiparticle spectra. 

2.1 Model 

To avoid encumbering the discussion with unessential details, we consider scalar 
QED, i.e. a charged scalar field <f> coupled to an external vector potential A^. 
Moreover, in order to simplify things even further, we neglect any kind of self- 
interactions among the scalar fields, and the coupling to the external electro- 
magnetic field enters only via the covariant derivatives, = 9 M — ieA^. Thus, 
the Lagrangian of this model is: 

C = (D^) (D"0)* - m 2 ^* . (1) 

In most of the considerations of this paper, we need not specify the precise form 
of the background potential A 1 *. Only in the final section, we work out the case 
of an explicit example of background electric field that leads to exact analytical 
results. 

2.2 Reduction formulas 

We assume that the initial state of the system does not contain any particles or 
antiparticles. However, because of the background field, transitions to populated 
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states are possible. Let us consider the following transition amplitudes, 

M m ,n({Pi}> (gJ) = ( Pi ' ' ' Pm out|O in >, (2) 

particles antiparticles 

from the vacuum to a populated state. The conservation of electrical charge 
implies that an equal number of particles and antiparticles must be produced, 
i.e. that this general amplitude is proportional to S mn - This transition amplitude 
can be obtained from the expectation value of time-ordered products of fields: 

/m n 
JJ d 4 x, (D xt + m 2 ) JJ &% (□„, + m 2 ) 

i=l j=l 

x (O out |T0(xi)---0(x m )0*(2/i)---0*(j/„)|O in ) . (3) 

Here the on-shell boundary condition is the vacuum one, i.e. — > \/p\ + m 2 
and — > \/ + m 2 for particles and antiparticles, respectively. This is ade- 
quate only if one chooses a gauge 2 in which the background field vanishes 
when time goes to +00. Because the conjugate <p*(yi) already takes care of an- 
tiparticle nature, qf should also be chosen to be positive. Note that, in principle, 
each field in this formula should be accompanied by a wave-function renormal- 
ization factor, Z^ 1 / 2 . However, since we do not include any self- interactions 
among the fields, these factors are equal to unity here and we can safely ignore 
them. 



2.3 Generating functional: definition 

All the physical quantities related to particle production in this model can be 
constructed from the squared amplitudes |A4 m , n | 2 . A very useful object that 
contains all this information in a compact form is the generating functional 
defined by [26, 28] 

00 ^ „ m n 

e ^/n d3 p«*(Pi)n d3 9j*(9j) 

m,n— i—1 j — 1 

(4) 

In this functional, z and z are two functions defined over the 1-particle momen- 
tum space (unlike what the notation may suggest, they are independent and 
not complex conjugates of each other). 

If one sets the functions z and z to constants equal to unity, one gets, 

00 

^[1,1]= P ^n, (5) 

m,n— 

2 It is always possible to find such a gauge if the electrical field vanishes when time goes to 
infinity, a necessary condition to be able to unambiguously define what we mean by "measuring 
a particle". If one insists on using a gauge in which is not zero when x° — » +00, one 
must replace the ordinary derivatives by covariant derivatives and the plane waves by gauge 
transformed plane waves in eq. (3). The mass-shell condition for and gP should also be 
altered by the non-zero background gauge field. 



M mt n{{Pi},{Qi\) 
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where P m ,n is the total probability to have m particles and n antiparticles in 
the final state. From unitarity, the sum of all these probabilities must be equal 
to one, hence 

^[1,1] = 1. (6) 

This is an important constraint on the generating functional !F[z,z\, that leads 
to significant simplification in the computation of inclusive observables. 

Assuming that this generating functional is known, one can obtain the single 
inclusive particle spectrum as 

dN+ _ 6F[z, z] 
d 3 p Sz(p) 

the single inclusive antiparticle spectrum as 

dN- _ 5T[z, z] 
d 3 q 5z(q) 

and the double inclusive particle-particle spectrum as 

dN++ <5 2 JP[z, z] 

d?p x d?p 2 5z(p l )5z{p 2 ) 

Other combinations of inclusive 2-particle spectra are given by 

dN 2 ~ _ 5 2 F[z, z] 
d 3 q 1 d 3 q 2 8z(q l )5z(q 2 ) 

Note that from their definitions, these two particle spectra are normalized so 
that their integrals over p and q are, respectively, 

where N^ 1 denote the number operator for particles and antiparticles in the 
final state respectively. In terms of the total probability introduced in eq. (5) 
we will see that (N+) = £ ro ,„ mP mn , (N~) = E ro ,„ nP mn , (JV+(JV+-1)> = 
S m n m { m ~^)Pmn, etc, for which one can find a justification in the appendix A. 

What these equations mean in the ++ and cases is that our 2-particle 

spectra are defined by summing over all possible pairs of distinct particles in 
every event. When summed over all momenta in a given event, this leads to 
N ± (N ± — 1) where N 1 * 1 is the multiplicity of particles (resp. antiparticles) in 
that event. Naturally, this requirement of taking distinct particles has no in- 
cidence on the H — case - since charged particles are always distinct from the 
corresponding antiparticles -, which explains the last of eqs. (11). 



(7) 

Z — Z — 1 



(8) 

2 = 2=1 



(9) 

Z — Z—1 



5 z T[z, z] 
d 3 pd 3 q Sz(p)6z(q) 



divr 
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2.4 Generating functional: computation 



Let us now proceed to the actual computation of the generating functional 
!F[z,z\. There is usually no closed form answer for this object. Since we are 
neglecting the self-interactions of the fields <f> in our model, however, this be- 
comes a much simpler calculation. It has been shown before [26, 28] that the 
generating functional T[z, z] is the sum of all the vacuum-vacuum 3 graphs in 
a slightly modified version of the Schwinger-Keldysh formalism [48], where the 
off-diagonal components G° + _ and G°_ + of the free propagator are altered by 
the functions z or z. Explicitly, the propagators read: 



As one can see, the off-diagonal free propagators are simply multiplied by z(p) 
and z(p) respectively. Given these propagators, the rules for calculating T[z, z] 
are straightforward: 

i. Draw all the vacuum-vacuum diagrams at the desired order. There are 
simply connected and multiply connected graphs. However, one can al- 
ways limit the calculation to the simply connected ones, and then expo- 
nentiate the result in order to obtain the full result that also includes the 
multiply connected ones. 

ii. For a given graph, sum over all the possible ways to assign + or — signs 
to the vertices. 

iii. A — vertex is the complex conjugate of a + vertex. Let us denote by 
V(A) the value of the coupling of 4>,4>* to the background field in a + 
vertex. Note that this 'potential' is not simply itself, since there are 
both a e(d ll (j))(f>* and a e 2 (fxj)* A^A^ couplings - however, we will not 
need its detailed expression in the following. The corresponding — ver- 
tex is V*(A) = —V(A) (this identity follows from the hermiticity of the 
Lagrangian) . 

iv. Connect these vertices with the propagators defined in eq. (12). 

Step i is trivial: there is only one topology of simply connected vacuum- 
vacuum graph in our model. These are the graphs made of a single closed loop, 
embedded with the background field A^, as illustrated in fig. 1. We must sum 
over the number of insertions of the background potential A^ (from zero to 
infinite insertions), and for each of these insertions we must sum over the type 
+ and — for the corresponding vertex. This double summation can be organized 
in blocks, as illustrated in figs. 2 and 3. 

3 Vacuum- vacuum graphs arc diagrams that have no external legs with respect to (j>. 
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p2 _ m 2 _ j £ 



G ( l_(p) = 27r9(-p a )z(p)S(p 2 -m 2 ) , 
G°_ + (p) = 2tt9(+p°) z(-p) 5{p 2 - m 2 ) . 



(12) 
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Figure 1 : Topology of the connected vacuum- vacuum diagrams that contribute 
to In T . The solid line denotes the free propagator of the charged scalar field <j> 
(the arrow indicates the direction of the flow of positive electric charge). The 
wavy line terminated by a circled cross denotes the background gauge potential 
V(A). Note that, in scalar QED, the background 'potential' is not simply 
itself, since there arc both a e(d l _ l (f>)(j)* A 11 and a e 2 (jxj>* A^A^ couplings. The 
wavy line represents the sum of these two contributions. 




Figure 2: Building blocks for the summation of field insertions having a fixed 
Schwinger-Keldysh vertex assignment. 




Figure 3: Block decomposition of the double summation over the number of 
background field insertions and the ± assignments at the vertices. 
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If we denote by T + the sum of graphs in the first line of fig. 2 and by T- the 
sum of graphs on the second line of the same figure, we can take the remaining 
steps ii, iii, and iv and the sum of the vacuum-vacuum diagrams contributing 
to In J 7 can be written as 



In T[z, z] = constant + 



1 

- tr 



T + G° + _T_G°_ + 



(13) 



where the trace symbol (tr (•••)) denotes an integration over the space-time 
coordinates (not represented explicitly in the formula) of all the vertices. The 
first term in this formula, that we simply denoted "constant" but did not write 
explicitly, is independent of z and z. It is made of all the graphs in which all the 
vertices are of type + or all of type — (and thus cannot contain z nor z since 
these come with the G± T propagators). In the second term of this formula, 
the index n represents the number of the block consisting of one H — and one 
— h transitions, and the factor 1/n is a symmetry factor since we can rotate the 
graph by one block without altering it. Note that the index n gives the order in 
z and z of the corresponding term 4 . 

It is trivial to perform explicitly the summation in eq. (13) to find, 



\nf[z, z] = constant - tr In 1 - T + G + _T_G 



(14) 



The constant term that we did not write explicitly can be determined without 
any calculation so that In 1] = 0, as required from unitarity (6). Therefore, 
we have 



F[z,z] = 



exp ^— tr In 


1 - T + G° + _T_G°_ + 


) 


exp ^— tr In 


1 - 


-T+G\_T_G a _ + 


Z — 


J 



(15) 



Although fairly formal, this formula contains all we need to know about the 
particle production by an external electromagnetic field in scalar QED. 

In order to simplify the subsequent discussion, let us restrict ourselves to 
background electric fields that do not depend on the position x. In this case 
it is always possible to choose a gauge in which the background vector poten- 
tial A^x) is also independent of x, and its Fourier transform is proportional 
to a delta function S(k) as far as its dependence on the spatial components 
of the momentum is concerned. In this particular case, the 2-point function 
T + G^_T-G°_ + has the same entering and outgoing momenta. In order to make 
more explicit the z and z dependence, let us introduce the notation: 



z(p)z(-p)L f 



(16) 



The important points here are that L does not contain z and z, and that the 
functions z and z carry the same momentum up to a relative sign 5 . With this 

4 We see explicitly here that the order in z of a given term is the same as its order in z, 
which reflects the fact that particles and antiparticles can be created only in pairs. 

5 Physically, the building block in eq. (16) is the amplitude squared for producing a single 
particlc-antiparticle pair. Since the background field is uniform in space, the total momentum 
of this pair must be zero. Hence the opposite sign for the momentum argument of z and z. 



9 



compact notation: 



cxp ( 


— tr In 


1 - 


zzL 


) 


exp 


(-tr In 


1 


- L 


^ 





From now on, it is simpler to calculate the trace in momentum space. Indeed, 
when the background potential is space independent, a unique momentum p 
runs around the loop. 

To close this subsection, let us mention a generic property of the trace that 
appears in eq. (17). Strictly speaking, when the background electric field is 
independent of the location in space, this trace exhibits a factor (2ir) 3 S(0) in 
momentum space. This factor should be interpreted as the volume V of the sys- 
tem 6 , and its presence is an indication that the particle spectra are proportional 
to the overall volume. 



2.5 Relation with wave propagation in the field 

So far, we have not attempted to calculate the object L that appears in the 
generating functional. Since it is built from the 7± , which are Feynman (time- 
ordered) propagators amputated of their external legs, it is clear that L is re- 
lated to the propagation of small fluctuations over the background electric field. 
However, knowing that L is related to the Feynman propagator is inconvenient 
for practical calculations because this propagator obeys complicated boundary 
conditions. In practice, one should try to rewrite L in terms of propagators that 
obey simpler boundary conditions, like the retarded propagator. 

Let us start from the equation that defines 7+ to rewrite it into the re- 
tarded quantities. The resummation that leads to T + can be summarized by 
the following Lippmann-Schwinger equation: 

r+ = v + vG° ++ r+ , (is) 

where V is the sum of the two couplings to the background field (the derivative 
coupling to a single and the non-derivative coupling to A^A^). We do not 
need to specify more what V is. Note that we could have written the equation 
in a slightly different form: 

T+ = V + T+G\ + V . (19) 

(This just amounts to starting the expansion from the other end-point of the 
propagator.) Concerning 71, it is sufficient to note that it is the complex 
conjugate of T + . 

6 In order to check this, one can make the background electric field slightly space dependent, 
so that it has a compact support in space. One sees now that all the integrals are finite, and 
that the single particle spectrum is proportional to the size of the region where the background 
field is non-zero. 
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The amputated retarded propagator T R is defined from the same equation, 
but the free Feynman propagator G\ + is replaced by the free retarded propa- 
gator: 

t r = v + vg r t r = v + t r g° r v , (20) 

where the free retarded propagator G° R (p) is defined as 

p z — m z + ip u e 

In order to express 7+ in terms of T R , the first step is to relate the Feynman and 
the retarded propagators. This is done via the following well-known relationship: 

G° ++ =G° R +p., (22) 

where p_ is a 2-point function whose definition in momentum space is 

p ± { P ) = 2ir6{±p°)6{p 2 -m 2 ). (23) 

(Note that is nothing but G\_ with z = 1.) From the above equations, we 
arrive trivially at 

(1 - VG° R - V P -)T + = (1 - VG° R )T R , (24) 

and subsequently at 

(l-(l-VG° R )- 1 Vp-)T + =T R . (25) 



Thus, we have 
Similarly, one can prove, 



T + = (l-T R p_)- 1 T R . (26) 



T + =T R (l-p_T R )- 1 . (27) 



In fact it is easy to confirm that eqs. (26) and (27) are equivalent by expanding 
the inverse quantity in terms oiT R p- in cq. (26) and P-T R in eq. (27). Taking 
the complex conjugate of eq. (26), we get, 

r_ = r; = (i - r^p-r'r; . (28) 

(p_ is purely real.) Multiplying eq. (27) by p_ on the right, we finally obtain, 

r +P _ = r R (i - p-rj- 1 p_ = r R P _ (i - r^r 1 . (29) 

Combining everything, we obtain the following expression for zzL: 

zzL - T+zp_T_z P+ = r R zp_{\ - T R p_)-\l - T*p_)- 1 T*zp + . (30) 
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Thus we have managed to replace all the Fcynman propagators in L by retarded 
ones. The price to pay for this transformation is that we have now an expression 
that is no longer bilinear in the propagators, but has terms at any order > 2. 
As we shall see now, this apparent complication actually disappears thanks to 
an identity reminiscent of the optical theorem. 

The Lippmann-Schwinger equation for T R is given in eq. (20). For T* , it 
reads: 



T* = -V - VG°*T* = -V - T*G°*V 



(31) 



Here, we have used the fact that, in a unitary theory (like scalar QED with a 
real background potential), we have V* = —V. Adding up the equations for T R 
and T* we first obtain 



t r +t: = -t:g o :v+vg o r t r 



(32) 



Finally, we can eliminate V in the right hand side of this equation, in favor of 
T R or T*. This leads easily to: 



T +T* 



Gt+G 



P+ 



(33) 



Note that this relation is a variant of the optical theorem 7 applied to a 2-point 
function. The left hand side is equal to the discontinuity of the 2-point function 
across the real energy axis, and the right hand side gives the expression of this 
discontinuity in terms of cut graphs. Thanks to this relationship, it is now 
straightforward to check that 



(i - r; P _)(i - 7>_) = 1 + T*p + T R p_ 



(34) 



By combining eqs. (30) and (34), we easily arrive at the following simplifi- 
cation; 



1 — zzL = 1 — T R zp_{\ + T* P+ T R p_)- l T*zp + 



= (i + r> + r R p_)- 1 

which leads to the generating functional: 



l-{zz-l)T R p_T*p + 



T\z, z] = exp ( - tr In [l - (zz - l)T R p_T*p + ^ ) . 



(35) 



(36) 



Here again, thanks to the fact that the background field is uniform, all the 
factors inside the logarithm share a single spatial momentum p. The function 
z has argument p and the z is evaluated at —p. 

At this point, we see that all the properties of the distribution of produced 
particles are determined by a single quantity, namely the amputated retarded 
propagator T R of a scalar particle on top of the background field A^. Before 



7 This is why the relation V* = —V, that is the manifestation of unitarity in this calculation, 
is crucial in order to obtain eq. (33). 
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going further, it may be useful to rewrite the combination zzT R p_T*p + with 
all the momentum dependence: 

zzT H p_r; P+ ] M = z(-q)p + (q) J z{k) T R (p,k)p.(k)T*(k,q) . (37) 

This formula is completely general. Note that T* is the same quantity as T R , 
in which all the retarded propagators are replaced by advanced ones and V(A) 
is replaced by V*(A). 

In eq. (37), the momenta k and q are forced to be on the in- vacuum mass- 
shell, since they appear inside the distribution p±. However, in the case of k, it 
turns out to be simpler to have a momentum variable (which we denote here by 
k) that obeys the mass-shell condition imposed by the non-zero gauge potential 
at x° — > — oo. Let us denote G^ the retarded propagator evaluated in the 
presence of the (constant) background field 

A^= lim A^x). (38) 

Since there is no electrical field at x° — > — oo, the gauge field is a pure gauge in 
this limit: 

A™ = d^ X (x) , (39) 

and the propagator is simply obtained by a gauge transformation from the 
vacuum propagator G° R : 

G™ {x, y) = e iex ^ G° R {x, y) e ~ iex ^ . (40) 

Let us now rewrite the combination T R p_T* that appears in eq. (37) in terms 
of the corresponding expressions 7^° and p°? which are naturally functions of 
the modified mass shell momentum k. This can be done by writing eq. (37) as 

T^T^T^KG^G-iGlr^p- ((Glr'G™)* ((G~)-iG R T R y . 

(41) 

Note that because 

(D™+m 2 )G™(x,y)=5(x-y), (42) 

where is the covariant derivative constructed with the asymptotic field A^ 
at y° — > — oo, and p°?{x,y) is still translationally invariant, it now projects 
momenta to the mass shell in presence of the background field at x° — > — oo. 
Since this gauge field is a pure gauge, it is easy to write the corresponding 
mass-shell conditions imposed by p+(k) and p°?(k): 

(k±eA°°) 2 = m 2 , (43) 

where the signs ± are to be chosen for the positive and negative energy solutions 
respectively. In a gauge where A = 0, as we shall chose later on, these read: 



Too p oo (T oo } » 



fc° = and fc° = -E% where = V \k + eA 00 ) 2 + m 2 . (44) 
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Therefore, eq. (37) can be rewritten as 

zT R p_T*p + ] = z(-q)p + (q) [ ^L Z (k)T™(p,~k)p™(k)(T™(k,q)y 
= z(p)z(-q)p + (q) J |r-(p,-fc)| 2 . (45) 

In the second line, we have changed k — > — k, and we have exploited the fact 
that T R x (p, —k) is proportional to S(p + k) in a uniform background field. Note 
that thanks to the constraints provided by the p+(q) and p+{k) factors, the 
energies k and q° are both positive. However, they obey different mass-shell 
conditions. The outgoing particle 8 energy q° follows the in-vacuum dispersion 
relation, while k° obeys the dispersion relation in the presence of the background 
field A°°. 

For practical calculations of T R yD (p 7 —k), it is best to relate this quantity to 
the Fourier coefficients of a plane wave propagating on top of the background 
field. Since is obtained by amputating the retarded propagator G R with 
(G° R ) _1 on the right and with (G^°) _1 on the left, we can immediately write: 

T™(P, -k) = J d A x jv* (U x + m 2 ) j d A y (D™ 2 + m 2 ) G R (x, y) 



lim [ d 



3 xe [ P- x (d X0 - iE° p ut ) Vk (x) . (46) 



In these formulas E° ut is the vacuum on-shell energy E° ut = \/p 2 + m 2 . Since 
the propagator G R (x,y) is a Green's function of the operator D 2 + m 2 (now, 
the covariant derivative is defined with the full background field, not just 
its asymptotic value in the past), 

[D 2 x +m 2 ]G R (x,y) = S 4 (x-y), (47) 

it is easy to check that rj k (x) obeys the following equation of motion: 

[D 2 x +m 2 ] Vk (x) -0, (48) 

provided that k obeys the negative energy mass-shell condition (43). In order 
to find the boundary condition when x° — > — oo for rjk(x), we can replace the 
full propagator G R (x,y) in eq. (46) by the propagator G°^{x,y) that resums 
only the asymptotic field 

Vk (x) = f d 4 y e iL y (D™ 2 + m 2 ) G™(x,y) 

e lhx . (49) 



8 Recall that we will need the trace of j^zT^p-T^p+J and thus q will be equal to the 
momentum of the produced particle p. 
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In order to obtain the final formula, we have used: 



(D™ 2 + m 2 )G™(x,y)=6(x-y) . (50) 

Thus, we see that the initial condition for r/k (x) is a plane wave, with a momen- 
tum k that obeys the mass-shell condition of eq. (43). 

In a uniform background field, we can simplify a bit the notations by writing 

T~(p, -k) = -2\E™\2k)H{ P + k) P P , (51) 

so that 

z*T~p~(T~)V+| =2E™ t (27r)H(p- q )z(p)z(-q)p + ( q ) \p p \ 2 . (52) 

The only quantity that we need to determine in order to fully solve the problem 
is the coefficient /3 p . This is obtained by solving the equation of motion (48), 
with a plane wave initial condition when x° — > — oo. We note that the initial 
plane wave is chosen as an antiparticle-like one in eq. (49) and projected into 
a particle- like one in eq. (46), the intuitive meaning of which will be clear in 
discussions in sec. 3. 



2.6 Multiparticle spectra 

Let us now use eqs. (37) and (52) in order to obtain results about the spectra 
of the produced particles. From now on, let us simply denote E° ut as E p in this 
section, for we have chosen the definition of (3 P so that E™ will never appear in 
the expressions. We shall wait for the next section where the difference between 
^m.out an( j physic interpretation of j3 p will be more explicit. The single 
inclusive particle spectrum is obtained as the first derivative of the generating 
functional with respect to z(p). We obtain 



z,z—l 



= (2^ l^l 2 • < 53 > 

It should be mentioned that, in accord with the definition (4), the functional 
differentiations with respect to z(p) and z(q) are not accompanied by (2ir) 3 . In 
the second, we have made the trace explicit. The final S(p— k) comes from the 
differentiation with respect to z(p). In the last line, we have performed the d 4 fc 
integration explicitly, and we have interpreted the (infinite) factor (27r) 3 i5(0) as 
the volume V of the system. Naturally, the spectrum of antiparticles is identical. 
For later reference, it will be useful to introduce more compact notations as 
follows: 

f P ^ { ^n P =\p P r 
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Note that f p has the interpretation of the occupation number for the produced 
particles of momentum p. This is clear if we integrate eq. (54) over the momen- 
tum and write it as 

<~ + > = /w^ (55 » 

where the properly normalized phase space measure d 3 pd 3 £e/(27r) 3 , or restoring 
the Planck constant d 3 pd 3 x/h 3 , is explicit. 

Let us now turn to the 2-particle spectra. For two particles, we obtain 

diV 2 ++ _ dN^_ dN+ 
d 3 p 1 d 3 p 2 d 3 p x d 3 p 2 
S 2 



Sz(Pi)Sz(p 2 ) 



tr [zzp+ T R p_ (T R ) zzp+ T R p_ (T R ) 



Z — Z — 1 



= S(Pi -P2)n P J Pl ■ (56) 

The uncorrelated part of the 2-particle spectrum shows up naturally in this 
calculation, and we have absorbed it in the left hand side. The right hand side 
represents the correlated component of the 2-particle spectrum. As one can see, 
particles are correlated only if they have identical momenta. By integrating the 
previous equation over p l and p 2 , and by using the first of eqs. (11), we obtain: 

(N+(N+ - 1)) - (N+f = J d 3 p n p f p , (57) 

or equivalently 

(AT+7V+) - (AT+) 2 = / d 3 p n p (l + f p ) = J f p (l + f p ) . (58) 

The form with an explicit integral over x is the one that would be applied to 
a system with a phase space density that depends (slowly) on the coordinate. 
The first term in the right hand side (the 1 in 1 + f p ) is the answer one would 
obtain for a Poisson distribution. That is, if the probability distribution is given 
by 

c-W(AO m 



m 



(59) 



which defines the Poisson distribution , one would have (AT + AT + ) — (N + ) 2 = 
(AT + ) (i.e. the variance and the mean are identical). Thus, the deviations 
from a Poisson distribution are contained in the term proportional to n p f p . 
Equation (58) indicates that these correlations are Bose-Einstein correlations, 
i.e. due to stimulated emission of particles in a single quantum state. 
For one particle and one antiparticle, we get: 

d7V 2 +" dN+ diVf _ S 2 



tr [zzp + T™p™{T™y 



d 3 pd 3 q d 3 p d 3 q 5z{p)5z(q) 

+ tr [zzp + T R x >p°°(T R °°yzzp + T~p~(T~)*] 
= 5(p + q)n p (l + f p ) . (60) 
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In this case, the correlation can only exist if the particle and antiparticle have 
opposite spatial momenta. Again, the integrated form of this equation reads 

(N+N-) - (N+)(N-) = J d 3 p n p (l + f p ) . (61) 

One may have wondered why the particle and antiparticle have correlations. 
This can be understood by the fact that, as we noted repeatedly, the pair 
of a particle and an antiparticle is created at once so that the particle and 
antiparticle production preserves the momentum conservation as well as the 
charge conservation. Therefore, the correlation (60) reflects the particle-particle 
correlation (56) with p x = p and p 2 = —q, which explains the delta function of 
spatial momenta in eq. (60). 



2.7 Nature of the distribution 

The results of the previous subsection suggest that the distribution of produced 
particles obeys the following properties: 

i. Two particles are correlated only if they have identical momenta. 

ii. A particle and an antiparticle are correlated only if they have opposite 
momenta. 

iii. In a given momentum mode, the distribution of produced particles follows 
a Bose-Einstcin distribution. 

Let us now present a more general justification of these results. Because 
the background electric field is uniform, a unique momentum p runs around the 
loop. By using results obtained in the previous subsections, we can rewrite the 
generating functional explicitly as follows: 

F\z, z] = exp (-V J In [l - (z(k)z(-k) - l) f k ] ) . (62) 

In writing eq. (62), as already discussed previously in this paper, we have as- 
sumed that the system is placed in a finite volume V. This is indeed necessary 
in order to have a finite particle production rate in a constant (in space) external 
field. A consistent quantization of the system in a finite volume requires one to 
specify boundary conditions 9 at the edges of V, which leads to the momentum k 
being a discrete variable. The continuum in k is recovered in the limit V — > oo. 
Switching now to a notation which makes this explicit, and remembering that 
d 3 fc/(2-7r) 3 = J^k we can wr ite the generating functional as 

^M=n uf \ . f > ( 63 ) 

j~ 1 + Ik - k k Ik 



Periodic ones being the most convenient. 
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where we use a compact notation; z& = z(k) and zp. = z(— fc). From this 
formula, one sees immediately that the distributions of produced particles in the 
various modes are totally uncorrelated, since the generating functional factorizes 
as a product of generating functions for single modes: 

F[z,z] = T\ Fk(zk,Zk) , Fk(z k ,z k ) = ——? — = — j~ ■ ( 64 ) 

x fe x 1 + fk~ z k z k fk 

By Taylor expanding this formula around Zk, z k = 0, it is easy to obtain the 
probability of having particles and antiparticles in the mode i, 

P( W ,K}) = IlT?f {ytj-J ■ ( 65 ) 

Given the occupation numbers fk (obtained from (3k by solving the equation 
of motion of a plane wave over the background field), this formula completely 
specifies the distribution of produced particles and antiparticles. Distinct modes 
are not correlated. In each mode, there must be an equal number of particles 
and antiparticles. The distribution of the particle multiplicity in the mode 
A; is a Bose-Einstein distribution of occupation number fk- A Bosc-Einstein 
distribution is in sharp contrast to a Poisson distribution (which would be the 
result in a complete absence of correlations), since its decrease at large rrik is 
much slower because of the absence of Wfc! in a Poisson distribution (see the 
denominator of eq. (59)). As a result, final states with many particles in the 
same momentum mode are more likely. 



3 Bogoliubov transformation interpretation 

We can interpret the results obtained in the LSZ derivation as a Bogoliubov 
transformation. To do this explicitly it is useful to switch to canonical quan- 
tization, which we shall review here shortly as the following manipulations are 
very standard ones. 

From the Lagrange density of eq. (1) one obtains the Hamiltonian of the 
theory (apart from the gauge part), 

H = J d 3 x nnt + ieA n$ - ieA ltf& + (m 2 - e 2 ^)$$t + (£>$) • (Z)$) t 

(66) 

where $ and n are operators in the Heisenberg picture satisfying the equal-time 
commutation relation, 

[*(*), II(y)] =iS(x-y). (67) 

It will be convenient for the following discussion to choose a gauge where Ao = 0, 
because in this gauge it is possible to directly associate the time dependence 
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of the wave function with the physical energy of the particle 10 . We shall thus 
work with the following Hamiltonian, 



H.js 



X 



nnt + m 2 $& + (£>$) • 



(68) 



The whole dynamics of the matter fields is determined by the equations of 
motion 



d $ = i[H, $] =n f 

d U = i[H,H] = (D 2 -m 2 )& 



(69) 



These can then be expressed as an equation of motion for $ only, but involving 
second order time derivatives. The important thing to realize is that because 
we are looking at a theory without self-interactions and coupled to a classical 
background field 11 , the equations of motion of the field operators are linear; 
they are in fact the same as the classical equations of motion for the fields. The 
solution to the retarded field equations therefore contains all the information 
about the relation between the field operators $ and II at a; — > — oo and 
x Q — ► oo. In the Heisenberg picture, knowing the relations between the field 
operators is equivalent to knowing the dynamics of the theory; in particular the 
whole probability distribution of the produced particles. 

We can now introduce the familiar decomposition of the field operators in 
terms of creation and annihilation operators. One can perform this decompo- 
sition in different bases of operators; in particular the ones that correspond to 
particles at x° — > — oo (the "in" states) or a; — > oo (the "out" states). Since 
these are just decompositions of the same operator 3> in different bases, one gets 
the equality, 



d 3 fe 

d 3 fc 

(2^)3 



■2E£ 



2E U 



out, A; 



(70) 



Here we have denoted the dispersion relations of particles and antiparticles in 
the in-state as E 1 ^ and that in the out-state as ££ ut . We assume that the 
background electric field is turned off adiabatically x° — > ±oo, which means 
that Ai approaches a constant value A™' ou . The dispersion relation for par- 

(fe + e A in ' out ) 2 and the one for antiparticles 



tides is then El' 



V 



I Consider for example a particle at rest in the vacuum: in the gauge = its wave- 
function is e — ,mx ° . Performing a time-dependent gauge transformation with the function 
Mx°/e will generate a (constant) gauge potential Aq = —M/e and change the wave-function 
to e — M ) x , which for M > m will seemingly look like a negative energy one. 

II Coupling to a quantum field would induce effective self-interactions through loop correc- 
tions. 
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£,in,out, _ ^/ m 2 _|_ Jfc _ e _A_in,out)2^ j n W riting eq. (70) we have used the sym- 
metry E k = E~^ k to write everything in terms of the particle dispersion relation 

£in,out = £ ,in,out,+ _ ^ ^ momcntum labcl fc rcfcrs to thc « canon i ca l" 

momenta, which are the variables describing the oscillation of the wavefunction 
in space. The momentum that is actually measured in a detector is the "ki- 
netic" one, which in this case is k + e A ln ' out for particles and k — eA ln ' out for 
antiparticles. We are keeping the notations rather general in this section. In the 
physical situation we are interested in, the only particles that are measured are 
the "out" -ones. A convenient gauge choice, and the one adopted in sec. 4 is then 
to take A° ut = 0, so that one need not distinguish between the canonical and 
kincmatical momenta for particles in the final state-it is enough to remember 
that E% is different from E% nt . 

The choice of basis in the decomposition cq. (70) of the field operator is 
determined by the boundary conditions for the functions <fj^ k (x) and (j)f ut k (x). 
When we require that they approach plane waves at asymptotic times: 





g— iE l £x°+ik-x 


for 


x°- 


-> — oo 


(71) 




— gi-Efc n x +ife-a! 


for 


x°- 


-> — oo 


(72) 


'out,fe( X ) 


— c -iEZ ut x°+ik-x 


for 


x°- 


-> +oo 


(73) 


out,fc( X ) 


= e i-E° ut :r +ifc-x 


for 


x°- 


-> +oo 


(74) 



the corresponding operators ai n ,out, ^in.out annihilate the in-state and out-state 
particles and antiparticles, respectively. Note that for further convenience our 
notation has been chosen such that the coordinate dependence in both (j> k (x) 
and 4> k {x) is e + and thus the usual negative energy plane wave e corre- 
sponds to (f)Z k (x). The canonical commutation relation for $ and II is satisfied 
if 12 

[ain,k,al ] = [6 in ,fe, b\ ] = [a out;fe , aj ut>p ] = [6 ou t,fe, &Lt,p] = (27r) 3 <5(fc - p) . 

(75) 

All the space-time dependence of the field operator $ is in the coefficient func- 
tions (f> in out k (x); the creation and annihilation operators are time- independent. 
Because the equation of motion (69) for $ is linear, the coefficient functions 
^in out k ( x ) mus t each independently satisfy the same equation. In fact the so- 
lution to the equation of motion rj k {x) introduced in eq. (48) is nothing but 

The relation between the field operators at x° — > — oo and at x° — > +oo is 
encoded in thc Bogoliubov coefficients. They are in the transformation matrix 
between the in- and out-basis functions. The solution for fe (x) is again a 
superposition of plane waves at x° — > +oo. If the background field depends 
only on time, thc modes of different k do not mix and we can introduce the 
Bogoliubov coefficients as the coefficients of this plane wave decomposition by 

12 Note that in here these commutation relations do not include a factor 2E^ as is conven- 
tional. The normalization used here is simpler in the case where E 1 ^ differs from E^ ut . 
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writing 

j™ x <*(*) = Jj£ («* e-"^- "^- + ^• "•) • (76) 



By noticing that fe (x°, — a;)]* satisfies both the same initial condition as 
(/>r k (x) and the same equation of motion, one finds the solution at x° — > oo 
that starts as a negative energy wave as 



J™J^ k W = \j^ (a* k e iE * Mx ° +ik - x +P k e-^" +ik -) . (77) 

This can also been seen using c~ lfc :E </>^ fe (x) = [c^' 3 ^;" fc (x)]*. We can now 
deduce the relation, 



<*(*) = \/^t («fc0+ t , fc (x) +^^- t>fc (x)J , (78) 



/ E m / \ 
^,*(*) = y ("feCut.feW +&# utifc (a0) • (79) 

In the general case of a space dependent background field the Bogoliubov 
coefficients arc not diagonal in momentum space. However, one can diagonal- 
ize the transformation matrix from the in- to the out-states, and our following 
discussion will equally well apply to the eigenstates of the more general transfor- 
mation instead of individual momentum modes. Inserting the decompositions 
(78) and (79) into eq. (70) one gets 

aout.fc = a k a in , k + Pk b\ n _ k , 
. ' . ' (80) 

KnUk= a -k b kk+P-k a in,-k ■ 

Consistency with the commutation relations (75) gives the normalization con- 
dition, 

\a k \ 2 - \[3 k \ 2 = 1 . (81) 

This normalization condition is a consequence of the charge conservation sym- 
metry of our Lagrangian 13 . The above relations can be inverted to give 

a m ,k = a* k a ut,k - Pk b f out _ k , (82) 
b in ,k = a*_ k b out , k - P-k aLt -k ■ ( 83 ) 



3 One can verify that eq. (81) implies 

Qout = Qin > wherc Q ~ e J (T^j- ( a l a fe ~ b k b k) 
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The vacuum is defined as the state that is annihilated by all the destruction 
operators (we shall drop the in- and out-labels for a moment), 

o fc |0) =6 fc |0) =0. (84) 

A properly normalized state with n particles of momentum k can be constructed 

as 

\n h ) = (-±=) |0>. (85) 



Particles of momentum k are counted with the particle number operator 

dN _ a k a k 



d 3 fe (2ir) 3 



(86) 



For example, the expectation value of the number operator on a state with one 
particle of momentum p is 

After setting up these conventions let us return to the problem of particle 
production. We consider the situation where there are no particles in at x° — > 
— oo, then the system is in the incoming vacuum state defined by 

ain,p|0 in ) = 6in,p|0 in ) = . (88) 

We are working in the Hciscnberg picture where there is no time evolution 
in the states, so the system stays in the |0i n ) state. But at late time a; — ► 
+oo particles are described by the "out" annihilation operators which do not 
necessarily annihilate |0; n ). In order to count the number of outgoing particles 
we must count the number of out-particles contained in the state |0i n ). For 
this it is useful to derive an expression for |0; n ) in terms of the "out" quantities 
flout.p, b out , p and |0 ou t)- That is, we have to solve, 

a in,p|0i n ) = {(X p dout,p ~ /^p^out,-p) l^in) = , 
&in,p|0in) = (a*L p 6out,p - /3-p%ut,-p) l°in) = • 

This above equation can be easily solved using the following ansatz; 

|0 in ) = C]Jc W [\ k al utk bl ut _ k )\0 out ) , (90) 
k 

where C is a normalization constant. Applying the condition (89) to the 
ansatz (90) we find the condition, 



n=i 



" E -fih («L, fe )"(Ct,- fe )" +1 ] lOout) = . (91) 

n=0 
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Shifting the summation variable by one in the first term gives 

00 in 

E [ X * a ** V ~ (*U k ) n (bU- k ) n+1 lOout) - , (92) 



n 

n=0 



which is satisfied by 



X p = y- 1 % . (93) 



a 



Now we can fix the normalization constant from 



(Oout I exp| A* a out , p 6 ut,-p } expj A p aJ ut p 6* ut _ p 1 1 ou t> 



1 , ,2 (M) 

= 1 + \Pp\ 2 , 



i-v 2 |a p | 



where we used the normalization condition (81). We then have the expression 
for the initial state in the following form: 



io in ) = na + i/3fci 2 r 1/2 ex P 



v "out.fc °out,-fe 

a fe 



|0out> • (95) 



This has a clear physical interpretation that the initial vacuum is a superposition 
of states with out-state pairs of particles with k and antiparticles with — k. 

Armed with the explicit expression (95) it is now straightforward to cal- 
culate, for example, single and double inclusive spectra by taking expectation 
values of the number operator. The spectrum of particles is 



dN+ 



^. <w ^«| ( y.^ l(1 + l / yVE-(^)"-^rl*r 

(96) 

This expression exactly coincides with eq. (53). The two particle spectrum 
(equal sign) is likewise 

diV 2 f+ f a lut,p 1 a out,p 1 a lut,p 2 a out,p 2 ^ut.p^out.P! \ . . 

\ U in| /0 so , ro 0{p 1 - p 2 ) , r~ |Uin) 



d 3 Pid 3 p 2 x 1 ^ (2tt) 3 (2tt) 3 ^ l (2tt) 

= (27r) 6 ^ in l a ° ut ^Pi a out,p 2 a out,p i aout,p 2 |0in) 

Note that we are explicitly subtracting the delta function contribution to agree 
with the definition (11). The variance is 



(N+N+) - (N+)(N+) = J d 3 p n p (l + /p) , 



(98) 
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which should remind the reader of the particle number fluctuations in a Bosc- 
Einstein system. It is actually easy to see that this is indeed what we have by 
looking directly at the probability distribution. 

In order to see more clearly the structure of the probability distribution we 
note that we can decompose the Fock space into a direct product of the Fock 
spaces for different momenta k. We then write an arbitrary state as a tensor 
product 

l*>=®l*>fc- (99) 

k 

It is convenient to group together particles of momentum k and antiparticles 
of momentum k under the same label k. Thus we define the vacuum states 
|0, 0)^ of the subspaces k as 

|0,0)fc = 6out,-fc|0,0) fc = 0. (100) 

We can now write |0 out ) as a tensor product of the vacuum states of the different 
modes, 

|0out) =®|0,0> fc . (101) 

k 

A state with m particles of momentum k and n antiparticles of momentum — k 
at x° — > oo is then denoted by 

""•">- fe) (w) |0 *- (102) 

By taking tensor products of the states \m, n)k for different k one can construct 
a complete basis for the Fock space. 

Applying this decomposition to eq. (95) we write 

|0i„) = (g) { (1 + \Pk\ 2 T 1/2 cxp \ V -^al nt k bl tt _ J |0, 0) fc j , (103) 
i. I L a k J J 



and expanding the exponential we get 

"Ifc 

_ U '^i i +■ U "nut U 



|0 in )=(g){(l + |/3 fc | 2 )- 1 / 2 £ ±(v- 1 £aU>bU-*Y''\W) 



m h =0 



(104) 



TOfc 



(1 + \Pk\ 2 )- 1/2 ]T (%) \m k ,m k ) k \ . (105) 



k K m h =0 x fc ' ) 

Thus we see that explicitly the "in" vacuum for a momentum mode k is a 
superposition of outgoing particle-antiparticle pair states. The amplitude for 
being in a state with m k pairs is 

M fe = fc (™ fc ,m fe |P fc |0 in ) = -=== (t±) , (106) 

V 1 + \0k\ 2 \ a kJ 
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where we have to introduce ¥ k , the projection operator to the subspace fc, to 
take the inner product 14 . The corresponding probability to have mj. pairs in 
the mode fc and any number of particles in the other momentum modes is 



P{m k ) = k(m k ,m k \F k \0i, 





( w \ 




f A ) 


1 + 1/3* I 2 




-i + A 





(107) 



This law for the probabilities characterizes the Bose-Einstein, or geometrical, 
distribution. Since the momentum modes are independent, we can then write 
the combined probability distribution as 



n*^> -n its (rfe 



(108) 



which is exactly the form at which we already arrived in eq. (65). 



4 Exactly solvable example 

So far our discussions and formulas are given in a rather general way. In what 
follows, we consider an example which is exactly solvable and demonstrate how 
the formulas in the preceding sections lead to the concrete evaluation of spec- 
trum of produced particle under an external electric field. 

The LSZ reduction method and the interpretation as a Bogoliubov transfor- 
mation both require that the asymptotic states at x° — > ±oo are well defined. 
This means that the external fields should be adiabatically vanishing so that we 
can define |0j n ) and |0 out ) without ambiguity. Thus, the imposed electric fields 
must be time-dependent, beginning with zero at x a = — oo, growing finite with 
increasing x° , and diminishing to zero again as x° — > +oo. 



4.1 Choice of the gauge potential 

It is known that the Klein-Gordon or Dirac equation under the following time- 
dependent electric field E = (0, 0, E(x )) is exactly solvable [3, 34, 47]; 

E(x°) = ^ . (109) 

[cosh(cjx )] 

This electric field exponentially goes to zero for \x°\ 3> In the limit of 

lo — > the electric field becomes homogeneous in time. We can choose a gauge 

14 The states |0j n ) and |"ifc , »rifc)fc live in different Hilbert spaces; one in the whole Fock 
space and the other one in its subspace fc. In order to take an inner product one therefore 
has to project out the fc component of the state |0i n ). Physically this projection means that 
we arc not measuring the other momentum modes of the state |0i n ) than fc. Thus eq. (106) 
gives the amplitude to have pairs in the mode fc and any number of particles in the other 
momentum states. 
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in which Aq = and the vector potential associated with the electric field 
(109) is A = -(0, 0, A 3 = Jdx°E{x°)) lb . We can immediately recover eq. (109) 
from E = — V^4° — 8qA. After the integration we find the Sauter-type gauge 
potential, 

r*° F 
A 3 {x) = / dy a E(y°) = - [tanh(wz ) - l] , (110) 

where we have chosen the integration constant so as to make A 3 (x) — > at 
x — > +oo. As we have discussed in sec. 3 this is a very natural choice in 
the present case. A constant ^3 amounts to a shift in the third component of 
the momentum p 3 , which is interpreted as a different frame choice. It is most 
natural to sit in the frame in which the particles and antiparticles measured at 
x° = +00 are at rest if their p 3 is zero. The gauge potential and the electric 
field are sketched in fig. 4. 




Figure 4: Sketch of the chosen gauge potential and the associated electric field 
as a function of time x°. The electric field has a peak at x° = whose height is 
E and width is specified by 1/ui. The gauge potential has an offset from zero by 
—2E/uj in the in- vacuum at x° = —00, meaning that the origin of p 3 is shifted 
by this offset. 



4.2 Solving the equation of motion 

An explicit solution for the Sauter-type gauge potential is already known, and 
so we will simply explain the necessary notation and then jump into the known 
expression of the solution. To make this paper as self-contained as possible 
we supplement the derivation in appendix B in more detail. Introducing the 
following notation, 

the equation of motion in scalar QED is given by the gauged Klein-Gordon 
equation as 

do - (d 3 - iAw[tanh(wx°) - l]) 2 - d\ + m 2 | <j>(x) = , (112) 



3 We work with the (+, — , — , — ) metric convention. 
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which is immediately read from the Lagrangian density (1). Because the equa- 
tion of motion does not have explicit dependence on x, we can factorize the wave 
function <j)(x) into the time-dependent part and the spatial plane-wave part, i.e. 

<Kx) = e ik x . (113) 

The time dependence is governed by the differential equation, 

do + (fc 3 + Aw[tanh(wx°) - l] f + k\ + to 2 ] ip k (x°) = . (114) 
We here change the variable x° into £, defined by 

£ = i [tanh(wa; ) + l] . (115) 

With this variable we can eliminate the hyperbolic function from the equation 
and express it only in terms of meromorphic functions. Moreover, £ has a 
convenient asymptotic behavior. We will later make use of 



1 for x° — > +oo 



1 + c 2wx ° J I e 2iaulx for x° 



-oo 



(116) 



and also 



+oo 

1 + c - 2uJX " J { 1 for x v -» -oo 



(l-O ta =(T^=wr) , Z!o^ + !! (117) 



to infer the plane-wave boundary conditions. 

Also, for concise notation we introduce the dimcnsionlcss energies in the 
same way as in Dunne's review [3]; 

fin pout 

where, with the time-dependent gauge potential, the energies in the in- and 
out-states are, respectively, 

(E%) 2 = (k 3 - 2Aco) 2 + k\ + to 2 , {E° k ut ) 2 = kl + k\+m 2 . (119) 

We note that the energy of antiparticles is given by E™ k in the in-vacuum and 
E'—k m l ne out-vacuum, respectively. Also we should explain our convention of 
the electric charge e. Our choice is as follows; particles are negatively charged 
and antiparticles are positively charged for e > 0. These are reminiscent of 
electrons and positrons in real QED. Therefore, in the above, noticing that 
k 3 = — fc 3 we see that the particle dispersion relation E k n starts with a larger 
longitudinal momentum than E^, which is understood as the deceleration by 
the Lorentz force in the direction anti-parallel to the external electric field. 
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Now we are ready to express the solution of the differential equation. One 
can express two linearly independent solutions in terms of the hypergeometric 
functions as 

^ 1) (0=r iAt (i-0- i ^i(l-i(A'+M+^),l+i(A'-M-^);i-2iM;0 ; (12Q) 

=^(l-0" i!y 2i ;, i(|-i(A'-M+^),|+i(A'+/.-^); l + 2i//; , 
where A' = yj\ 2 — 1/4 with A defined in cq. (111). 



4.3 Asymptotic behavior of the solution 

Now we need to take an appropriate linear combination of the two solutions 
written in cq. (120), so that the the boundary condition like eq. (48) can be 
fulfilled. In fact, it will turn out that these solutions already satisfy a simple 
plane- wave boundary condition. To make this explicit we should use eqs. (116) 
and (117) together with the general property of the hypergeometric function; 
2-F\(a, 6; c; £ — > 0) — > 1. Therefore, in the limit of x° — > — oo (i.e. £ — > 0), we 
see, 

(£ ~* °) - e- 2i ^° = c-^" , Vi 2) K -> 0) -> c 2i ^° - c 1 ^" . (121) 
In accord with the identification defined by eqs. (71) and (72) we have 

Cfc^e*"". ^.fc^e*", (122) 

The quantities necessary to compute the particle and antiparticlc production 
are obtained as the coefficient of these solutions in the x° — ► +oo limit-when 
decomposed in terms of the cjy^ nt k . To this end, it is necessary to know the 
limit of 2-Fi (a, b; c; £ — > 1). One must be, however, very careful when taking 
this limit because a, 6, and c are complex numbers in this case. Thus, it is 
convenient to make a transformation in the argument from £ to 1 — £, which is 
possible by means of the following mathematical identity; 

2 F 1 (a,b; c; x) = ^ 2 Ji(o,6; 1 - c + a + 6 ; 1 - x) 

r(c-a)r( c -6) ^ 123 ^ 

+ r(C ^^p ) (1 - sr—^itc - a, c - b ; l+c-a-6 ; 1 - x) , 
which leads to an alternative expression of the solution, 

Vi 1) (0=r i ' i (l-0" i,/ ^2F 1 (i-i(A'+M+^) ) i+i(A'-M-^);l-2i^; 1-0 
+ r if(i_^)iM fc2 F 1 (i+i(A , -Ai+i/),J-i(A , +A i -i/);l+2ii/; 1-0, (124) 

and 

4 2) (e)=r(l-e)- i ^ 2 F 1 (i-i(A'-/i+ I ,),i+i(A'+ A1 - ; ,); l-2hv; 1-0 

+ ^(l_^ Bfc2 J? 1 (l +i ( A ' + Ai + 1/ ) } l_ i ( A '_^_ 1/ ). l + 2il/ ; 1-0, (125) 
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where we defined, 

T(l - 2i / u)r(-2hv) 



A 



k 



T(i - i(A' + ,x + u))T(i + i(A' — n — v)) 

d* = r(i-2i M )r(2i^) 1 J 

fc r(i + i(A'-M + ^))r(i-i(A' + M -^)) ' 

Here we note that, if we take a naive limit of x — > 1 using a textbook formula 
on 2-Fi(q, b ; c ; 1), it would miss the second term in the identity (123) because of 
the factor (1 — x) c ~ a ~ h which is zero if a, 6, c are real but is an oscillatory finite 
function if a, b, c are complex. At this point it is straightforward to deduce the 
asymptotic plane- wave forms at x° — > oo, that is, £ — > 1; 

- 1) - A fe c- 2i — ° + B* k e 2i ^ ° = A fe e- iE * tx ° + B* k e iE °^^ 



(127) 

</4 2) K -> 1) - ^ e 2i — ° + B fc e- 2ii ^° = A* e i£ * ut * + B fc e"^ . 



Comparing the above behavior with the Bogoliubov transformations (76) and 
(77), we obtain the Bogoliubov coefficients as 



/ 7^0llt / 7^0Ut 

^V^F 4 *' pk = yW Bk ' (128) 

4.4 Particle spectrum 

Before addressing the concrete expressions for the particle spectrum, we list all 
the necessary formulas to proceed with the calculations. The Gamma function 
generally satisfies, 

T(l + z) = zT(z), r(l - z)T(z) = , (129) 

sin(7rz ) 

from which we can derive the following useful relations, 

|r(ia)| 2 = . - - , |r(i + kO| 2 = — - , |r(i + ia) 



asinh(7ra) ' sinh(7ra) ' 2 cosh(7ra) 

(130) 

Then, after some algebra, we reach, 

, ,2 _ cosh[7r(A / +£+ v)\ cosh[7r(A' - fj, - v)\ 

l(lkl ~ sinh(27r/i) sinh(27r^) - { ' 



and 



2 cosh[7r(A' — /x + z/)] cosh[-7r(A' + \x — v)] 



' /3fc ' sinh(27r/x) sinh(27w) ' 

We note here that, using cosh(a + b) = cosh(a) cosh(fe) + sinh(a) sinh(6) twice, 
we can easily check that 

|« fc | 2 - \p k \ 2 = 1 , (133) 
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which is consistent with the condition (81). Now that we have \(3k\ 2 explicitly, 
we can get the general probability distribution which is characterized only in 
terms of |/3fe| 2 . The single inclusive spectrum, for example, is 

dN+ _ V cosh[7r(A / - p, p + v p )\ cosh[7r(A / + \i p - v p )\ 
d 3 p ~ (2tt) 3 sinh(27r/ip) sinh(27w p ) ' ^ ' 

From this expression we can get the occupation number f p which is obtained by 
removing the volume factor V/(2tt) 3 of the single particle spectrum. Once f p is 
given, the whole probability distribution is known as discussed in the previous 
sections. We plot f p as a function of p 3 in the unit of m± = \Jp\ + m2 m n g- 5. 
In drawing fig. 5 we fixed mj_, and set the electric field to the value E = Trrn^/e- 
which is sufficiently strong to create particles in view of the standard expression 
of the Schwinger mechanism-, and then we vary the time scale uj. 




Figure 5: Time scale dependence of the produced particle distribution as a 
function of p 3 . As u> — * the distribution approaches eq. (135) and f p extends 
between p 3 ~ —2eE/u> and p 3 ~ Oas seen from the curve for u — tot in the 
figure. In contrast, with increasing uj, the result approaches eq. (136) which 
spreads wider than the small-w case with the distribution center located at 



-eE/uj. 



Now let us consider two extreme cases. First, we take the constant field limit 
(u) — ► 0), which make the above expression as simple as follows; 



V 



diV+ 
d 3 p ~* (2tt) 



■ exp 



tt(p 2 ± + m 2 ) 
AeE 



1 



1 



(u 0) , (135) 



where p = cup 3 / (2eE) taking a value in the range of —2eE/uj < p 3 < (i.e. 
— 1 < p < 1). In the outside region, p 3 > or p 3 < —2eE/uj, the result is zero 
in the to — * limit. 



30 



This distribution of particles in the p 3 direction corresponds to the momen- 
tum distribution as a result of the deceleration by the electric field while it is 
imposed. If ui is small, the electric field lives long, and the produced particles 
are pushed down by the Lorentz force along the p 3 direction for longer time. 
Note that this range bounded from —2eE/ui to zero coincides with the range 
of A 3 changing between x a = ±00. One might expect that the total number 
of produced particles diverges in the case of a constant (in time) electric field. 
This is indeed true; the result of the p 3 integration is nearly proportional to the 
integration range 2eE/u when u) is small enough, which is divergent as 1/w-i.e. 
as the time during which the external electric field is non-zero. 

Next, we shall take a look at the opposite limit, i.e. a short-pulse limit, 
u — > 00. In this limit eq. (134) is reduced to, 



To have non-zero value the electric field eE should be larger than u>E° ut . Even 
though there is no exponential suppression 16 , as compared to the small-w case, 
the resulting f p is significantly suppressed by large u> for a fixed maximal 
strength of the electric field, as is apparent in fig. 5. 

Let us finish this subsection with a comment on a related work by Cooper and 
Nayak [49]. In this paper, the authors study the Schwinger mechanism for the 
pair production of charged scalars in the presence of an arbitrary time-dependent 
background electric field. In particular, they calculate the pair production rate 
(i.e. the number of pairs created per unit of time) 17 , and conclude that "the 
result has the same functional dependence on E as the constant electric field E 
result with the replacement: E — ► E(t)" . However, the comparison of the two 
limits in eqs. (135) and (136) indicates that the time-dependent case is unlikely 
to be given by a mere replacement E — ► E{t) in the time-independent result. 

5 Conclusions 

We have calculated the multiplicity of produced particles under a spatially con- 
stant but time-dependent electric field in scalar QED using a formalism based 
on the LSZ reduction formula. We defined and computed the generating func- 
tional of the particle and antiparticlc distribution, from which we determined 
the whole distribution of the production probability. We found that particle 

16 When the electric field has a fast time dependence, the perturbative process 7 — ► <jxj>* can 
produce particles. 

17 One may question whether this is a well defined concept, since the presence of a time 
dependent external field makes the definition of proper particle states ambiguous. One may 
argue that the only quantities that can be defined unambiguously are those where measure- 
ments are done only after the external field has died out. From the point of view of the 
Bogoliubov transformation, the transformation coefficients are uniquely determined from the 
asymptotic plane-wave time-dependence. Of course one may use some working definition of 
the Bogoliubov coefficients at arbitrary intermediate time, but such a treatment implicitly 
assumes a quasi-static approximation. 





(136) 
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production in one momentum mode follows a Bose-Einstcin law, reflecting the 
statistics of scalar particles. We have also derived the same results by means of 
a Bogoliubov transformation on creation and annihilation operators. 

A natural question may arise as follows; what is the distribution not in 
spinor (ordinary) QED instead of scalar QED? Our discussions in this paper 
were quite simple because we focused only on spin-zero scalar particles. In the 
case of spinor particles we need to deal with the spin structure, which brings 
unessential complication in, though the generalization is straightforward. We 
shall here just mention that the distribution of fermionic particles in the same 
momentum and spin mode follows a Fermi-Dirac distribution. We can confirm 
this immediately by replacing the commutation relation by the anticommutation 
one in the derivation based on the Bogoliubov transformation. The transforma- 
tion (80) is unchanged since this originates from the asymptotic behavior of the 
wavefunctions. The normalization condition (81) should be |afc| 2 + \ flk\ 2 = 1 
then to preserve the anticommutation relation of the transformed operators. It 
is easy to show that \f3 p \ 2 gives the occupation number f p , and using the above 
normalization condition we can arrive at the probability distribution, 

P({m s ,k}) = - f s ,k) ( JZT7~ ) ' < 137 ) 

where s refers to the spin and m Sj fc takes the values or 1. Equation (65) (or 
cquivalcntly (108)) and the above (137) are our central results. We see that, 
if we drop higher orders than the quadratic terms in fk for fk <C 1, all of 
the Bose-Einstcin, Fermi-Dirac, and Poisson distributions are trivially reduced 
to identical answer; I — fk for no-particle production and fk for one-particle 
production. The difference emerges at the quadratic order, and it is notable 
that the difference remains no matter how small fk is. For instance, the two- 
particle production probability in the same mode is /|/2 if the distribution is a 
Poisson one, /| if a Bose-Einstein one, and zero if a Fermi-Dirac one. In other 
words one must properly take account of the quantum statistical nature if the 
multiparticle correlations are concerned. 

In the final section of this paper, we have revisited a known exactly solvable 
example of time-dependent electric fields. The time-dependence of the Sauter- 
type potential is actually ideal to think of the Schwinger mechanism; since we 
can unambiguously define the asymptotic states in the infinite past and future. 
This property of adiabatically vanishing external fields is necessary to make the 
discussion of particle production meaningful. 

Using the exact solution we took two extreme limits of constant (lo — > 0) 
and short-pulse (lu — > oo) electric fields. In the constant case we found that 
f p distributes almost uniformly over the range —2eE/uj < p 3 < 0. Therefore 
the total (integrated) number of produced particles diverges as 1/ui that is 
interpreted as the time duration for which the external electric field is imposed. 
In the short-pulse case, on the other hand, f p is a double-peak structure and the 
minimum in-between is located at p 3 = —eE/u. In practice the latter would be 
useful because it is more difficult to sustain larger eE/u in the laboratory. 
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A Multiplicity distribution 



In some cases, one is interested only in the overall number of particles and 
antiparticles in the final state, but not in their distribution in momentum space. 
For such observables, one can define a simpler generating function that does not 
contain any information relative to the momentum of the produced particles: 

= E S£ / fl d3 Piil d3 1j \M m ,n({P i }AQ i })\ 2 ■ (138) 

m,n— i— 1 7 — 1 



Note that this is also equal to: 

G[u,u] 



E - 

m,n— 



(139) 



where P m ,n is the probability to have exactly m particles and n antiparticles 
in the final state. Obviously this new generating function can be obtained 
from the generating functional ZF\z, z] by setting the functions z(p) and z(p) to 
constants respectively equal to u and u: 



G[u, u] = T[z(p) = u, z(p) = u] 



(140) 



Thanks to this relationship, one can obtain quantities such as those defined in 
cqs. (11) as ordinary derivatives of For instance, 



/ 



j3 J3 dN 2 + 



dG[u, u] 



d 3 p 1 d 3 p 2 



du 2 



u,u—l 



= J2 ro(m-l)P m ,„. (141) 



m,n— 



The last equality is obtained from eq. (139), and is the justification for the right 
hand side in the first of eqs. (11). 



B Detailed derivation of the solution 

With the variable change from a; to £, simple algebraic procedures lead to 
d a = 2wf(l - 0% and 8% = 4w 2 £(l - £)[£(! - £)d 2 + (1 - 2£)<%], from which 
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we can rewrite the equation of motion (114) in the following form; 

£(1 - Off + (1 - 2£)d c + M 2 r 1 + ^ 2 (1 - 0" 1 - A '] MO = • (142) 

In what follows we will explain how to solve this differential equation. Before 
finding the analytical solution, from this form of the equation we can already 
confirm that the asymptotic behavior of the particle solution is e^ 1 ^ x at 
x° — > — oo and e ±lE ^ x at i° ^ +oo as it should be. We shall pick up 
the most singular terms out of the differential equation (142), which gives 

[^| + ^ + M 2 r 1 ]4 n (0 = 0, (143) 

around £ = (i.e. x° — > — oo). It is easy to find the solution of this equation as 
^ n (0 = £±v ~ e ± iE k n x . In the same way we can extract the behavior around 
£ — ► 1 (i.e. x° — ► +oo) from the singular terms; 

[(i - oa| - a e + i, 2 (i - O" 1 ]^ (0 = o - ( 144 ) 

leading to Vfe Ut (0 = (1 - ±iu - e^ iE h* x ° . 

Now let us return to solving eq. (142). Because we have seen the boundary 
condition, it is convenient to factorize the plane-wave pieces as follows; 

#(0=^(1-0-^(0, (145) 
then we can find the equation that tpk(£) should satisfy as 

{£(1 - 0a| + [1 - 2i/i - (-2i M - 2w + 2)Z]d e 

(146) 

- (-i/i - w - W + l/2)(-i/x - \v + i\' + l/2)|^ fe (0 = . 
Here we recall that the hypergeometric differential equation, 

|a;(l - x)dl + [c - (a + b + l)x] d x - ab}f{x) = , (147) 
has two independent solutions given by 

fV>(x)= 2 F 1 (a,b; c;x), / (2) (x) = x l ~ c 2 F l {a+l-c, b+l-c ; 2-c; x) . (148) 
Therefore, by the identification of 

a = i - i(A' + // + v) , b = i + i(A' - fj, - v) , c = 1 - 2i/i , (149) 
we finally arrive at the solution (120). 
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